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Abstract. - Using a two-fluid model for viscoelastic polymer solutions, we study analytically 
fluid transport driven by a transverse, small amplitude traveling wave propagation. The pumping 
flow far from the waving boundary is shown to be strongly wave number and viscosity depen- 
dent, in contrast to a viscous Newtonian fluid. We find the two qualitatively different regimes: 
In one regime relevant to small wave numbers, the fluidic transport is almost the same as the 
Newtonian case, and uniform viscoelastic constitutive equations provide a good approximation. 
In the other regime, the pumping is substantially decreased because of the gel-like character. The 
boundary separating these two regimes is clarified. Our results suggest possible needs of two-fiuid 
descriptions for the transport and locomotion in biological fiuids with cilia and flagella. 
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J>~>' Introduction. — Active fluid transport or mixing in 
a complex fluid finds many examples on scales of cellular 

i__ibiology [IH^. In tlie fluid transport in respiratory or di- 
gestive organs [J] , large arrays of beating cilia on a surface 
Organize metachronal waves, in wliicli neighboring cilia 

,^^beat in a constant phase lag, to pump mucus fluids [5HI]. 
During development of an embryo, flows generated by ar- 
r,ays of rotating nodal cilia dictate the left-right symmetry 

T— I breaking in the placement of organs [5]. In plant cells, 
cytoplasmic streaming, which is responsible for rapid in- 
tracellular transport and mixing, is driven by the motor 
protein myosin that moves along filamentary actin at the 
cell periphery . 
^ . Hydrodynamic thrust generated by surface distortions 

kj>ife also used to propel many prokaryotic and eukaryotic 
cells [TU]- A sperm cell beats its single flagellum and 
^ims through mucus in the female mammalian repro- 
ductive tract [TT|. Microorganisms such as Paramecium 
or Volvox colonies possess ciliated surfaces, and propel 
themselves with metachronal coordination. As an example 
of locomotion in contact with a substrate, amoeboid cells 
crawl accompanying the protoplasmic streaming driven by 
periodic contractile cell body deformations [12]. 

All these examples involve small-scale, slow flows to 
which the zero Reynolds number hydrodynamics is ap- 
plied. In addition, fluid media can also be viscoelastic, 
such as a gel-forming character of mucus and protoplas- 
mic sol in amoeboid cells. While the basic principles of 
swimming at zero Reynolds numbers in a Newtonian fluid 



have been established for years [T H51IT5HT5] , the role of vis- 
coelasticity on flagellar motility or transport is still elusive. 

The importance of nonlinear effects of a fluid elas- 
ticity has been addressed in recent studies. Employ- 
ing various constitutive equations with fading memories, 
Lauga investigated a sheet with a given wave pattern and 
showed that the fluid elasticity always reduces its swim- 
ming speed [inillT]. Fu et al. extended the findings of 
Lauga to a beating filament, and discussed the changes in 
the beating pattern and of the swimming direction |18j . 
Teran et al. studied the problem numerically and demon- 
strated the enhancement of the swimming ability for large 
tail undulations [T^. All these works employed uniform 
constitutive equations for non-Newtonian fluids such as 
Stokes-Oldroyd-B models. However, in constrast to a 
Newtonian fluid acting as a uniform viscous background, 
viscoelastic media, such as polymer solutions and gels, are 
essentially dynamically strucutred materials with charac- 
teristic time and length scales. For example, transient 
network structures are speculated to be responsible for 
a mysterious viscosity-dependent swimming of some spi- 
ral bacteria pDH24j . To our knowledge, no theoretical at- 
tempt has been made to combine a mesoscopic description 
for complex fluids and a theory of active fluid transport 
and locomotion. 

In this letter, employing a two-fluid description for vis- 
coelastic fluids [5SHlHli we consider the two-dimensional 
hydrodynamic flow fleld driven by a planar flexible bound- 
ary with a transverse propagating wave of small ampli- 
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Fig. 1: Schematics of a waving sheet immersed in a viscoelastic 
fluid. The sheet (boundary) is infinitely extended and is fixed 
in the laboratory frame (except its undulating motion). The 
system is assumed to be uniform along j/-direction. A surface 
disturbance is a traveling wave of phase speed c and a net 
pumping fiow is generated at inifinity. 



tude |13j . and discuss the induced pumping flow at infin- 
ity (see fig. [IJ. This net flow corresponds to a swimming 
speed in the lab frame when the sheet is allowed to move 
translationally (TOlIIS]. We show that the hydrodynamic 
pumping flow generated far from the waving boundary 
is strongly wave number and viscosity dependent. For 
large wave numbers, the pumping transport is substan- 
tially decreased compared to the Newtonian case, and its 
asymptotic value is determined by the elastic properties 
of the medium. A mesoscale fluid property is crucial. The 
diffusive relaxation of the viscoelastic network stress and 
its frictional coupling to the viscous Newtonian solvent 
effectively acts as a sink of the velocity field, leading to 
the decrease of the far field pumping flow. Our analysis is 
thus complementary to the previous studies |16H18j , where 
spatially uniform viscoelastic stress relaxations were pos- 
tulated. 

Model. — In general, the hydrodynamic flow profile 
at the far field is rather insensitive to its near-filed de- 
tails. Therefore, many essential features and symmetries 
of flagellum propulsion or pumping with planar and heli- 
cal waves may be captured by the present minimal model. 
To proceed, we give a profile of the flexible boundary as a 
traveling wave of amplitude b 

z{x,t) = eh sm{kx ~ Lot) , (1) 

where we have introduced the small dimensionless expan- 
sion parameter e, which we will set unity eventually. The 
wave number k is taken positive without loss of generality. 
The displacement, eq. ([1]), is independent of a load, and 
the surrounding fluid satisfies the no-slip dounary condi- 
tions on the deforming sheet surface. 

We describe our model viscoelastic medium as a poly- 



mer network with mesh size coupled to a viscous New- 
tonian solvent with viscosity rj via a friction constant T, 
estimated as F ~ f?/^o PSVES] . The polymer network 
is described continuum elastic medium with shear 

and bulk Lame coefficients /i and A. Neglecting all iner- 
tial contributions unimportant to our aim, the linearized 
equations of motion for the network displacement u and 
the fluid velocity v are ^511^ 

AtV^u+(A + /^)V(V-u) = F(u-v), (2) 
-Vp + ryV^v = -F(u-v), (3) 
V- [(l-<?i)v + 0u] = 0, (4) 

where the dot represents the time derivative, and the in- 
compressibility of the whole fluid is required in eq. (|4|). 
Since the volume fraction of the polymer (p is sufficiently 
low even for entangled solutions, eq. Q is simplified as 
V • V = [551 [22]. We identify the elastic stress tensor 
of the network as ctP = /i[Vu -I- (Vu)^] -I- AV(V • u) and 
define a viscoelastic force f p = V • cr^ F(v — u). For an 
entangled polymer solution, the elastic network is actually 
formed only transiently; the elastic stress can be relaxed 
through disentanglements of polymer chains [^niSO] . This 
effect is most easily incorporated by adding a term r^^f^ 
to the dynamical equation for P derived from eq. ©-(HI): 

F = --P + /iV2v+^V2F + ^^V(V-P), (5) 
r F F 

where r is the characteristic stress relaxation time [?7l[5D] . 
Note that eq. ([5]) recovers the original two-fluid gel model 
for r — >■ oo [1SJ[15], while for F — >■ oo it reduces to a 
linear viscoelastic Maxwell model [3TJ[32]- Decomposing 
the vector P into its longitudinal and transverse parts 
as fP = il + where = V-2V(V • P), and = 
[1 — V~^VV] • P, we split eq. ([5]) into two independent 
equations. Transverse part of Stokes equation, ([3|), reads 
= V^v_L -|-f^. It is easy to see that a zero shear viscosity 
of this medium is given by 

■Hp = V + Tfi. (6) 

Perturbative analysis. — We introduce the stream 
functions, v = {d^ip, —9x4') and = (dzX, —dxX)i where 
dx denotes the partial derivative with respect to a;, so that 
they automatically satisfy the divergence-free condition. 
The longitudinal force, f^^, is by deflnition a gradient of 
a scalar potential: = — V0. The viscoelastic force is 
reconstructed from the scalar functions as = ~dx4> + 
dzX and /f = —dz(f> — dxX- The boundary conditions on 
the sheet surface read 

VV'(a;, e&sin(fca; — = ebco cos{kx — ujt)ex, (7) 

iP{x,ebsm{kx-ujt),t) = 0. (8) 

The first equation is the usual no-slip conditions for the 
velocity field v. The second one requires that the net- 
work and the solvent must move together on the surface 
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because no penetration of polymer chains across the sheet 
is allowed. 

Due to the linearity of the governing equations, we re- 
place time derivative with — and obtain the equations 
for the functions '0, x and 0: 



= TyV^VV + V^x, 



(9) 



= [v2-4H]v2x, o = [v2-K^H]v20, (10) 

where the inverse of dynamic screening length is given by 



1 + t'o — iujT 



1 — iujT 



i^A^)^ 72 ' «:||('^) = 72 ■ (11) 



The transverse and longitudinal viscoelastic lengths ap- 
pearing in eq. PT|) are 



and 



(A + 2/^)t 



r 



(12) 



where uq — T^iji] ~ rip/rj — 1 is the rescaled solution 
viscosity. Following the standard procedure, we expand 
■0, X and (f) in powers of e as t/i = eijji + e'^ip2 
X = ^Xi + ^^X2 + ■ ■ ■ and </> = e(f>i + e^02 + • • • • The 
governing equations at each order keep the same form as 
Eqs. © and pH)) . and the boundary conditions at 0(e) 
order are V?Ai = buj cos{kx)ex , f^i = —dx4>i + dzXi = 0, 
and f^i ~ ^dz4>i — dxXi = at z = 0. From the reflec- 
tion symmetry in the problem: £—>■—£ [3l[T3], there is no 
net flow and all viscoelastic force vanish at infinity. The 
first-order solutions take the following form: 



[Cie 



c.c, 



(13) 



(14) 



and 



^1 = [£:^e-V'='+«iM^ + (Fiz + Gi)e-'=^]e*('=^-'^*) + c.c, 

(15) 

where "c.c" stands for complex conjugate, and 
Ai, Bi,Ci, Di, El, Fi,Gi are complex numbers. Six 
relationships between the integration constants, two from 
eq. ([5]), and four from the boundary conditions, are found 
as Ai + ijk\{uj)Ei = 0, 2/irfcFi + (1 - iLxJT)Hi = 0, 
Ei + Gi = ibLj/(2k), y/k^ + Kl{Lu)Ei - Fi+ kGi = 0, 
k\{u)Ai + kHi + ikCi = 0, and 
kAi + kHi + iy^fc2 _^ Kjj(tj)Ci = 0, where Hi= Bi + iDi. 

Solving these algebraic equations uniquely determines ipi 
and thus the velocity field Vi. 

Next we proceed to determine the second order solu- 
tion. Looking at eqs. ([7]) and the boundary condi- 
tions at O(e^) read V'(/'2 = — 6sin(fcx — ujt)ydz^pi and 
= — 6sin(fca; — ujt)dz{f at z = 0. At this second or- 
der, a net flow at infinity is expected, which implies that 



time-independent terms should be added in the solutions 
for -02, X2 and 4>2- 



X2 



-Kj^(0)2 



= Woe 



-re||(0)z 



and 

V'2 



+ [£;2e-V4FT^iMz ^ ^^^^^ _^ Q^^^-2k.^^2^{kc,-ut) 

(18) 



c.c. 



where wcwcw^o, sq, k\{0) = (1 + i^o)/Ci ^'^'^ '^^(O) = 
(l-t-;/o)/C^ (sec eq. (|lip ) arc now real numbers. (The con- 
stants 742-G2 are complex numbers, but arc irrelevant to 
our aim below.) The solutions, (fT6)) - ((T8)) . satisfy dzX2 ~^ 
and 9z02 — ^ at z — >■ 00, which ensure that the polymer 
network and the solvent move together at infinity. 

The net pumping fiow u is the non-zero fiuid velocity at 
infinity: 

u ^ lim dziJ2 = e^so, (19) 

Thus the remaining task is to determine sq from the 
boundary conditions. After some calculations, we arrive at 
a uniform hydrodynamic flow at infinity of (setting e = 1) 



b^ujk 



[l + A'(fc,L.)], 



(20) 



where un = b'^ujk/2 = {bk)^c/2 is the result for Newto- 
nian fluids flrst obtained by Taylor [13]. The viscoelastic 
correction A' is the real part of the function A given by 

see eq. 

Equations ([20)1 and (PT|) constitute the main result in 
this letter. In flg. [21 u/u^ is plotted as a function of 
fc^x for varying wt, with the Poisson ratio of the network 
o- = set 0.2. 

Results and Discussion. — As seen in fig. [21 the 

pumping fiow is always smaller than that in a viscous 
Newtonian fluid. The most important feature is that it 
strongly depends both on {k,uj) and on the viscoelastic 
parameters. In particular, there are two qualitatively dif- 
ferent regimes. In small wave number or short viscoelastic 
length regime, the viscoelastic force "diffusion" in eq. ([5]), 
r^^/^V^fP, becomes negligible compared to r~^P, and the 
model behaves like a linear viscoelastic Maxwell model. 
Looking at flg. [H the pumping speed in this regime is un- 
changed from the viscous Newtonian case, i.e., u/uiq — > 1 
for all LOT, in agreement with the previous studies |31j . 
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A(fc,a;) = 



2k 



(21) 




ered as shown in the previous studies P^HTS] . 



100 



Fig. 2: The rescaled pumping speed, u/un, as a function of 
the rescaled wavenumber fc^ for different rescaled frequencies 
u}T=Q, 10,100, 1000. The rescaled viscosity and the Poisson 
ratio of the network are set vq = jir /-q = 100 and a = 0.2. 



In constrast, for large wave numbers, the stress propaga- 
tion is dominated by the shear elasticity, and the whole 
medium acts mainly like a cross-linked gel. The diffu- 
sive term, F^^/iV^P, now balances with /iV^v, leading 
to fP ^ Fv (thereby directly suggesting u ~ 0). The sol- 
vent flow through the virtually immobile rigid network is 
strongly damped via the friction (just like a flow through 
a porous medium [24]), and the far field hydrodynamic 
flow is accordingly decreased. Taking fc^x — > oo limit in 
eq. (j2ip . we obtain the asymptotic pumping speed 



Moo _ A + /,t(l + 2ri/rip) 



(23) 



This limit would mostly be relevant to dense viscoelastic 
cases, where ry/r^p <^ 1 and ^ Co- In such cases, we 
would obtain u/un ~ (A + n)/{X + 3/x), which indicates 
that the pumping fiow is ultimately determined by the 
elastic parameters only. 

Note that, while in fig. [5] we resaled the wave number k 
with the viscoelastic length the characteristic length 
separating these two regime is not simply but is gen- 
erally WT- dependent. In fig. [3l we show the "diagram" in 
the {k,uj) plane. Here the boundary line separating the 
two regimes is identified as the point at which u decreases 
the half amount of the full decrease in fig. [31 that is, where 
u = (l + 7ioo)/2. From this diagram, we see that increasing 
uj (for a given k) enlarges the uniform viscoelastic regime, 
where the Newtonian-like response, u/un ^ 1, is obtained 
within the linear analysis. This result for the first time 
clarifies in which condition the use of spatially uniform 
viscoelastic constitutive equations is justified. Note that 
in this regime, nonlinear elastic effects have to be consid- 
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Fig. 3: The diagram of the pumping response on the (fc, u) 
plane. Two qualitatively different regimes are found. The 
boundary line between the two regimes is determined as de- 
scribed in the main text, for different value of the solution 
viscosity: i^o = 10 (dotted, green), = 100 (broken, blue), 
and vq — 1000 (solid, red). 



With the scaling of the friction constant, F ~ v/^o, 
we obtain k^± ~ ^/vo{k^o)- The viscosity dependence 
for a fixed (fc,aj) is thus known from u/un = 1 + 
^'(v^(^Co); '^''")i which is plotted in fig. [Jj (assuming 
that r is independent of I'o). As a general tendency, the 
pumping speed descreases with ~ ilp/rj, in agreement 
with the previous theoretical and experimental observa- 
tions [Tni[IH]. For a very large solution viscosity, the vis- 
coelastic length, C i^y^'^S.o, also becomes very long, 
i.e., fc^x 3> 1, then the asymptotic value oi u/un is again 
given by cq. (|3^ . In stark contrast to the previous stud- 
ies [HJ , we predict the strong viscosity dependence of the 
pumping or swimming speed within the linearized vis- 
coelastic model. We stress that this is a direct consequence 
of the mesoscale structure generic to entangled polymer 
solutions. 

In the present two-dimensional geometry, the viscoelas- 
tic stress tensor can be given as — —0 + dxdzO, 
= -0 - dxd,0, and a^, = af^ = - ^(9^ - 01)9, where 
the function 9{x, z) satisfies = 2%. The total stress 
tensor is the sum of the viscous and the elastic contribu- 
tions: cr = — pl + ?7[Vv+(Vv)-^]+cr^'. Plugging the above 
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Fig. 4: The rescaled pumping speed, u/un, as a function of the 
rescaled viscosity, fo = rjp/ri — 1, for fc^o = 1 (with the relation 
k(,± = ^Juok^o), and for different ujt indicated on the figure. 



expressions for into cr, wc obtain 

(Jxx = -{p + 4>) + d^d,C, (24) 
(Tzz = ^{p + <t>)-d,d,C (25) 

a..=a,. = -\{dl-dl)(:. (26) 
where C = ^' + '^'ni' satisfies 

V^V^C = 0. (27) 

The work done by the waving surface on a surrounding 
viscoelastic fluid per unit time (the rate of working), W , 
is obtained by integrating the product of the velocity and 
the total force on the fluid over the surface of the sheet A: 

W = I v{cr-n)dA. (28) 

J A 

Transforming the surface integral in eq. (|28p into the vol- 
ume integral via the Gauss's theorem and using the me- 
chanical balance V • cr = (obtained by adding eqs. ([2]) 
and dni)), we flnd W = J{dkVi)aikdV . Plugging eqs. ((24)) - 
into this, and using Vx = dz^' and Vz = —dxip, we 
arrive at 

W = S + J dxj dz\7^ij\7^C, (29) 

where 5* is the surface integral given by 

S = I dx{dxi^){dx~dz){dzC)\z=o- (30) 

One can solve eq. (j27p if appropriate boundary conditions 
for the stress function (^{x, z) are given depending on the 
problem at hand. For example, for a freely swimming 
sheet, the total force acting on the sheet along the moving 
direction is zero: 

(j) ex ■ cr{x,z{x,t),t) ■ n = 0, (31) 



where fi is the unit normal to the surface, and the integral 
is over the surface of one period. Once (^{x, z) is known, 
the work W is calculated according to eqs. ([^ and (|5n)) . 
details of which will be published separately l33| . 

Finally, we briefly comment on the cross-linked gel (or 
solid network) limit, t — oo. This limit can be formally 
taken in eq. (j21[) without any mathematical singularity. 
The resulting function behaves quite similar to eq. (|2ip . 
with the characteristic length and time scale now replaced 
by (rj/ry^^ ~ ^0 and rj/^, respectively [29]. Note, how- 
ever, that cross-linked gels do not flow. In the gel limit, the 
no-slip conditions, eqs. ([7]) ([5]), are inapplicable. While the 
no-slip boundary conditions for the network displacement 
u require to drag the whole solid material along the pump- 
ing direction, it is impossible because of progressively in- 
creasing restoring forces. General slip boundary condi- 
tions, as well as convective nonlinearities, should be taken 
into account, as demonstrated for swimming in gels |34j . 

Summary. — By analyzing the viscoelastic two-fluid 
model, we have studied fluid transport driven by a flex- 
ible sheet with a given traveling wave deformation. We 
used the linearized hydrodynamic model, but appropri- 
ately included non-local nature of the viscoelastic stress 
relaxations. In contrast to a Newtonian fluid, the hy- 
drodynamic pumping flow generated far from the waving 
boundary is strongly wave number and viscosity depen- 
dent. In particular, the two qualitatively different regimes 
of the pumping response are clarified. In one regime, the 
Newtonian-like response is obtained within the linear anal- 
ysis, while in the other regime the pumping is substantially 
decreased due to the medium's gel-like character. The 
boundary line separating the two regimes is shown on the 
plane spanned by the actuation wave number k and fre- 
quency uj. The diagram illustrates when the spatially uni- 
form viscoelastic constitutive equations, such as a Stokes- 
Oldroyd-B model used in the previous studies [T5Hl9| . do 
provide good approximations. 

We conclude this paper with remarks on future exten- 
sions of this work. It is important to explore effects of 
nonlinear convective terms omitted in our viscoelastic con- 
stitutive equation. For finite amplitude undulations, this 
contribution may become substantial [TnillH]- More real- 
istic geometries should also be considered for a detailed 
comparison with cilia and flagella beatings. In particular, 
for free swimming, a crucial role of the undulating tail has 
recently been pointed out in ref. [121. Dynamical shape 
changes responding to the hydrodynamic forces should 
also be allowed. For those aims, numerical approaches 
would be more suitable. Effects of meso-strucures in a 
viscoelastic medium upon the active pumping could be ex- 
perimentally tested using artificial-cilia-mounted microflu- 
idic devices developed in refs. [551I5S] . Macro-scale experi- 
ments, similar to the one presented in ref. |37| . would also 
be informative. 

After we have submitted this paper for publication, we 
learned of the work by H. Fu, V. B. Shenoy, and T. R. 
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Powers, in which swimming in gels are studied using the 
two-fluid model [Ml. 
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